Water Chemisorption and Reconstruction of the MgO Surface 
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The observed reactivity of MgO with water is in apparent conflict with theoretical calculations 
which show that molecular dissociation does not occur on a perfect (001) surface. We have per- 
formed ab-initio total energy calculations which show that a chemisorption reaction involving a 
reconstruction to form a (111) hydroxyl surface is strongly preferred with AE = —90.2 kj mol -1 . 
We conclude that protonation stabilizes the otherwise unstable (111) surface and that this, not the 
bare (001), is the most stable surface of MgO under ambient conditions. 



PACS numbers: 82.65My 68.35Md 82.65D 71.10+x 



Magnesium oxide has long provided a prototype for 
the study of surface structure and chemical reactions of 
oxides. Naturally occurring MgO, known by its mineral 
name of periclasc, is not a common crustal mineral, but 
its simple structure makes it an excellent example for the 
investigation of mineral surface chemistry. 

Reactions at mineral surfaces are responsible for much 
of the chemical change which occurs in the Earth's crust. 
Weathering reactions control the erosion of rocks and the 
consequent evolution of surface topography thus provid- 
ing an opposing mechanism to the more dramatic process 
of mountain building. Aqueous reactions in sedimentary 
basins are responsible for the diagenetic processes which 
transform unconsolidated sediments into rocks. In this 
work we study the nature of a simple mineral surface 
when exposed to an aqueous environment and the chem- 
ical interaction of water with that surface. This is both a 
prerequisite to studying the interaction with aqueous so- 
lutions and a tractable first step towards ligand-exchangc 
reactions in more complex silicate minerals. 

We have performed experiments on single-crystals of 
MgO prepared with high-quality (001) faces which were 
reacted with acidic solutions. The experiments and re- 
sults are reported in detail elsewhereu, the main feature 
being the development of an altered surface layer. Elastic 
Recoil Detection Analysis (ERDA)cl shows protonation to 
a depth of 900 A with a H/Mg ratio close to 2 giving a 
probable chemical composition of magnesium hydroxide. 
Indeed brucite (the mineralogical name for Mg(OH)2) is 
the most common alteration product of periclase in the 
natural environment^ and well-crystallized intergrowths 
of brucite on periclase have been reported]. 

The initial stage in the reaction is hydroxylation of the 
surface. MgO has the cubic rocksalt structure with (001) 
cleavage planes. This is the moat stable surface and is 
the only one seen experimentallyO. The simplest possibil- 
ity for a hydroxylated surface is obtained by dissociating 



a water molecule and placing the OH group above each 
magnesium ion and the H above each oxygen of the (001) 
surface (see Fig. |l|a) as postulated by Coluccia et alu. 

Some striking hydroxylation experiments were repor- 
ted by Jones et al. who studied surface roughening 
on (001) faces of nanoccvstalline MgO in a transmis- 
sion electron microscopeu. The remarkable affinity of 
MgO for water is demonstrated by their in situ obser- 
vation of hydration-induced surface roughening over 10 
minutes under vacuum with Ph 2 o < 10~ 5 Pa. The 
presence of surface hydroxyl groups on MgO powders 
exposed to HaPi has been demonstrated by infra-red 
spectroscopyu'Bli3. Hydroxyls are clearly distinguishable 
from physisorbed molecular water by the HOH bend- 
ing mode which disappears above 100°C, while the OH 
stretching mode persists even above 500° C. Furthermore 
there is complete monolayer coverage of the surface by 
hydroxyls, as shown by microgravimetry measurements! 

Despite these observations, the most reliable the- 
oretical calculations predict that water molecules do 
not dissociate on the (001) surface. Scamehorn 
et a/Jlil calculated the energetics of the reaction 
{>MgO} + H 2 ^ {>Mg(OH) 2 }Bto form the hydroxy- 
lated (001) surface using periodic Hartree-Fock methods. 
They showed that dissociative chemisorption is energeti- 
cally unfavourable and that physisorption of intact water 
molecules is preferred. This was confirmed b«,a Car- 
Parrinello ab initio molecular-dynamics studyli3 which 
investigated the dynamics of a water molecule at a MgO 
(001) surface. No dissociation occurred. Experiments 
performed by Jones et alu also led to the conclusion that 
perfect (001) surface sites are not protonated. 

In summary, water demonstrably chemisorbs onto 
MgO but trustworthy calculations show that H2O mol- 
ecules should not dissociate on the only known stable 
surface. 

Several authors have proposed that water dissociates 
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FIG. 1. (a) Hypothetical hydroxylated (001) surface of MgO, (b) (0001) surface of Mg(OH) 2 , (c) our postulated (111) 
hydroxyl surface of MgO. The top three O-Mg-O layers of an oxygen-terminated MgO (111) surface have the same structure as 
Mg(OH)2 (0001). The (111) hydroxyl surface may be equivalently constructed by protonation of an oxygen-terminated surface 
or by hydroxylation of a magnesium-terminated surface. 



instead at low cooped inated sites such as steps, corners 
and other defectsBtj and dissociation-was demonstrated 
computationally at steps and corners! 1 2 tl 14 l. However one 
might expect this to simply saturate the defect sites. Dis- 
sociation at defects does not explain, the monolayer cov- 
erage of the surface by OH groupaJ, the whole-surface 
roughening observed in the TEM3 nor the progressive 
transformation of the entire (001) surface and incipient 
bulk hydroxylation observed by i.r. spectroscopy!'™ and 
in our dissolution experiments. 

In order to account for these phenomena we propose an 
alternative structure for the fully hydroxylated surface. 
The progressive formation of a bulk hydrated layer sug- 
gests consideration of the structural relationship between 
the oxide and the crystalline hydroxide. Mg(OH)2 is trig- 
onal and is composed of layers of Mg 2+ and OH - ions in 
(0001) planes (see Fig. |b). This "basal" plane is also the 
cleavage plane, yielding a stable type II hydroxide sur- 
face. The trigonal symmetry means that the basal plane 
is incommensurate with the cubic (001) MgO surface, but 
closely resembles the threefold symmetric (111). Fig. |l| 
illustrates that protonation of the (111) creates a surface 
with the same structure as the Mg(OH)2 (0001) cleavage 
plane. This strongly suggests that protonation may sta- 
bilize the MgO (111) surface, and that this may be the 
hydroxylated surface observed by i.r. spectroscopy. 

Supporting evidence is provided by jvell-crystallizcd 
natural growths of Mg(OH)2 on MgOd and from the 
inverse, reaction, the dehydration of Mg(OH)2 to form 
MgOE_3. In both cases an epitaxial relationship exists 
between the two phases such that the [0001] axis of the 
hydroxide is aligned with with the oxide [111]. 

We have investigated the stability of the (111) hydrosjd 
surface using total-energy pseudopotential calculations^. 
The density functional theory formulation of quantum 
mechanics is solved within the local density approx- 
imation— .(using the parameterization of Perdew and 
ZungertZI) by conjugate-gradient minimisation of the to- 
tal energy with respect to the valence electron wavefunc- 
tions. Optimized non-local pseudopotentials are used in 
the Kleinman-Bylander separable representation. The 
system of ions and electrons is subject to periodic bound- 



ary conditions which allows the use of a plane- wave basis 
set. The latter gives an accurate representation of the 
crystalline environment. Just as importantly, it permits 
analytic forces to be computed allowing for relaxation of 
the ions to their minimum energy configuration. This 
type of calculation has been used to study defect en- 
ergies hx MgOEa, OH groups as substitutional defects in 
quartzlij, reconstruction of the silicon fcU.1) surfacecJ and 
dissociation of CI2 at a silicon surfaced. 

One difficulty in using plane-wave pseudopotential 
methods for oxides has been that the tightly-bound oxy- 
gen 2p electrons require a high energy cutoff making the 
calculations expensive to perform. Recent advances have 
dramatically improved the convergence properties so that 
oxide calculations are now routineE3. In this work we 
used a new optimized oxygen pseudopotential which gives 
complete convergence of the electronic energy (to 0.05 eV 
from an energy of 909 eV) at a cutoff of 500 eV. Magne- 
sium was also represented by an optimized pseudopoten- 
tial and hydrogen by a pure coulombic potential. 

As a check on the accuracy to be expected we per- 
formed calculations on bulk MgO and an isolated water 
molecule, the initial reactants. The results agree closely 
with experimental values and are summarised in table Q. 



Property 



Calculated Experiment 



MgO 



H 2 



a (A) 
K (Mbar) 
K' 

bond length (A) 
HOH angle (°) 



4.217±0.001 
1.62±0.02 
4.19±0.08 
0.966 

103.9-104.1 a 



4.211 
1.603±0.003p 
4.15±0.10B 
0.9572 
104.52 



TABLE I. Calculated properties of bulk MgO and isolated 
water molecule: MgO lattice parameter ao, bulk modulus K 
and K' = dK/dP; H2O molecular bond length and angle. 
a The value depends on the size and shape of the supercell 
and the orientation of the molecule, indicating a very small 
interaction between a molecule and its periodic images. 

To represent a surface using periodic boundary condi- 
tions we performed calculations on a slab of MgO in a 
periodic cell of larger dimension leaving a vacuum sep- 
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arating periodic images of the slab. The surface energy 
is given by A_B surf = (-E s i a b - E hu i k )/a 2 . To test conver- 
gence with respect both to the number of layers in the 
slab and the vacuum space we calculated the energy of 
the (001) MgO surface using 4, 6 or 8 atomic layers in 
the slab and different c dimensions of the supercell leav- 
ing between 8 and 16 A of vacuum. The periodic cells 
had dimensions ao/v2 x ao/v2 x na Q ,n — 1,4,6. The 
lattice parameter ao was fixed at the experimental value 
of 4.2117 A. The k point set contained 8 points for the 
long cells and 32 or 48 points for the compact cell used 
for the bulk calculation. In every case relative energies 
were computed using equivalent sets to cancel basis-set 
size errors. The total energies were well-converged as a 
function of k-point sampling, to better than 0.05 eV out 
of 909 eV. The co-ordinates of the ions in the outer layers 
of the slab were relaxed to their minimum energy config- 
uration, while the inner 2 layers were constrained at their 
bulk separation. The surface Mg 2+ ions moved inwards 
by 1% and the O 2- ions by 0.1%.in agreement with pre- 
vious Hartree-Fock calculationsEi The resultant surface 
energies are all between 1.093 and 1.103 J m -2 which 
shows that a 4-layer slab with 8 A of vacuum gives ener- 
gies converged to a precision of better than 0.01 J m -2 . 
This result is compared with previous work in table |l| 



Technique 



AE suli (J 



DFT/LDA this work n 
DFT/LDA using Gaussian basis setEJ 
Periodic HF a ,-, 
Pair Potentials^ 
Experiment 13 



m J 
1.10 
1.16 
1.43 
1.07 
1.04-1.20 



TABLE II. Calculated MgO(OOl) surface energies and 
comparison with previous measurements and calculations. 
a These calculationsE] usedra-basis set containing only s and p 
orbitals. Birkenheuer et allEi showed that including d orbitals 
decreased their A_B suf -hfrom 1.32 to 1.16 J m -2 . 
Summarised by To^ andrefatherein . 

The test calculations establish that errors in the en- 
ergy differences due to cell-size and basis-set effects are 
less than 1 kJ mol -1 . This leaves one significant source of 
systematic error, the LDA, which is .known to underesti- 
mate molecular dissociation energies! 2 -! We estimated an 
upper bound for this error of 40 kJ mol -1 from a calcu- 
lation of the reaction MgO + H 2 # Mg(OH) 2 (brucite) 
whose energy is well known . Our main result is quite 
robust to an error of this magnitude. 

The chemisorption energy is fully determined given 
the structure and state of the reactants and products 
by AE = product - J2 ^reactants- The appropriate ini- 
tial states are the bare (001) surface of MgO and ice 
(since this is a zero-temperature calculation) and the fi- 
nal structures are the hydroxylated (001) and the (111) 
hydroxyl surfaces of Fig. |l|. Scamehorn et al. showed 
that the fully hydroxylated (00<L) surface has lower en- 
ergy than if partially hydratedO so we performed cal- 
culations for only that state. Calculations on the hy- 
droxylated (001) surface were performed in a supercell 



V^ao x v^ao x 4ao containing 4 layers of Mg and O ions 
plus 2 surface layers of hydrogen or hydroxyl. Simple 
electrostatic considerations indicate that full hydroxyla- 
tion must also be favoured for the (111) surface: this is 
the only means of achieving a non-polar surface. 

Ideally calculations on the reactants and products 
should use equivalent supercells to achieve cancellation 
of basis set errors in the computed AE. This rules out a 
direct calculation for the bare (001) ^ (HI) hydroxyl 
reaction as the symmetries differ and no cell can ac- 
commodate both structures. However the reaction may 
be split into two stages {>MgO} (001) ^ MgO bulk and 
Mg0 buik + jj 2 q ^ {>Mg(OH) 2 } (111) whose partial en- 
ergies sum to the desired result. Each stage can be com- 
puted using cells of the appropriate symmetry since bulk 
MgO is commensurate with both. The first stage is sim- 
ply the surface energy. The second used trigonal cells 
with a = ao/%/2, c = V3ao and 24 k-points for bulk MgO 
and a = ao/%/2, c = 2V3ao with 12 k-points and con- 
tained 5 layers of Mg and O ions in a hydroxylated slab. 
The outer layers were relaxed in all these calculations. 

The initial state of water was based on calculations 
of an isolated molecule in supercells equivalent to those 
used for the bulk MgO calculations but of twice the lin- 
ear dimensions to minimize interaction between periodic 
images. To this we added the experimental value £©£ the 
sublimation energy of ice, —47.35 ± 0.02 kJ mol -1 ! 2 -! 

The results are listed in table III. As expected the 



hydroxylation of the (001) surface is energetically un- 
favourable. Our hypothesis that hydroxylation stabi- 
lizes the (111) surface is confirmed, with a AE of - 
90.2 kJ mol -1 with respect to ice, or (by using the re- 
sults of Scamehorn et aliLH for the physisorption energy) 
-117.3 kJ mol -1 relative to physisorbed water. This pro- 
vides a ready explanation of the spectroscopic observa- 
tions of surface hydroxyls and the observed reactivity of 
MgO with water manifested as rapid roughening of the 
(001) surfaces of microcrystals. The chemisorption reac- 
tion must involve a reconstruction of the (001) surface. 



Initial 


Final 


AE (J m -2 ) 


AE (kJ mol -1 ) 


(001) 


(001) 


+0.78 


+41.6 


bulk 


(111) 


-0.59 


-31.3 


(001) 


(111) 


-1.69 


-90.2 



TABLE III. Computed energies for the chemisorption re- 
action of water with MgO, formally {>MgO} initial + H 2 



{>Mg(OH) 2 f 



The initial states are bulk MgO or the 



(001) surface plus ice at K. The final states are the hy- 
droxylated (001) or (111) hydroxyl surfaces. The (111) sur- 
face energies are expressed per unit of the original (001) sur- 
face area assuming a ratio of v3:l. A periodic HF study 
found a significantly more positive AE for the hydroxylation 
of (001) of 77-90 kJ mol -1 fil. The measured enthalpy of wa- 
ter chemisorption on MgO at 613 Pa and 543 K is in the range 

— 113 to —189 kJ mol -1 Lj. Approximate thermodynamic cor- 
rections to AE at K for the (111) hydroxyl give A_ffs43 ~ 

— 138 kJ mol -1 , consistent with those experiments. 

We conclude that the (111) hydroxyl, rather than the 
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bare (001) is the normal surface of MgO under ambi- 
ent environmental conditions. A (001) MgO surface will 
chemisorb water and reconstruct except under ultra-high 
vacuum or high temperature. Dehyekoxylation is exper- 
imentally observed only under UHVI13: indeed hydration 
reactions occur under lesser vacuum in the TEM at a par- 
tial pressure of water of < 10~ 5 Pa. Upon heating dehy- 
droxylation begins at 200° C but is gradual with residual 
hydroxyls persisting until over 700°CElll3. 

This solves a number of experimental puzzles. 1) The 
H 2 # D 2 exchange reactioa, is catalysed by MgO at tem- 
peratures as low as 78 Kc3. Structural surface protons 
provide the necessary exchange siteO. 2) It may also re- 
solve discrepancies between experimental and theoretical 
energies of adsorption of CO onto MgOcJ. The appropri- 
ate surface for molecular adsorption is the (111) hydroxyl 
surface, not the bare (001). 3) Vermilyea showed that the 
dissolution rates oil-MgO and Mg(OH)2 are identical over 
the pH range 2-5L3, an observation easily explained by 
the almost identical surface structures. 

Identification of the actual reconstruction pathway is 
beyond the scope of this paper, but a plausible mecha- 
nism must account for both the dissociation of the water 
molecule and the transport of surface Mg 2+ and O 2- 
ions. Previous ab initio calculations have shown that 
water molecules dissociate, at surface defects, particu- 
larly steps and cornerstSEJ, corroborating the observa- 
tion thatpthe reconstruction is more rapid if the surface is 
damagedlj. The activation barrier for ion transport may 
be estimated given that the inverse reconstruction, the 
annealing of {lOOUIacets on a (111) surface, is observed 
to occur at IOOOKEj. The equivalent thermal energy is 
12 kJ moU 1 , rather less than the chemical energies of 
the hydroxylation reaction. It is also possible that the 
barrier for Mg 2+ transport away from the defect site is 
reduced by hydration with either the product OH ion 
or an intact H2O molecule. This would expose low co- 
ordinated oxygen sites able to dissociate incoming water 
molecules and thereby continue the process. 

Stabilization of a type III polar surface by protonation 
or hydroxylation to form a non-polar type II surface is 
unlikely to be unique to MgO. We would also expect re- 
constructive chemisorption to occur in other fee binary 
metal oxides such as NiO. 

The consequences of the stability of the hydroxyl sur- 
face over the bulk are substantial but rather harder to 
predict. It does provide a driving force for the layer hy- 
dration observed by us and others which apparently leads 
to a topotactic oxide-hydroxide transformation. However 
a much better characterization of the hydroxylated layer 
is needed and further studies to establish its composition 
and structure are under way. Whatever its nature it is 
demonstrably a vital precursor stage in the dissolution of 
magnesium oxide and may also prove relevant to under- 
standing recrystallization and precipitation reactions. 
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